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Abstract 

In this article we describe the novel method to construct fundamen- 
tal solutions for operators with variable coefficients. That method was 
introduced in [11) to study the Tricomi-type equation. More precisely, 
the new integral operator is suggested which transforms the family of the 
fundamental solutions of the Cauchy problem for the equation with the 
constant coefficients to the fundamental solutions for the operators with 
variable coefficients. 
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1 Introduction 

In this article we describe the novel method to construct fundamental solutions 
for operators with variable coefficients. We also give a brief survey of some re- 
sults obtained by that method, which was introduced in [TP to solve the Cauchy 
problem for the Tricomi-type equation. Later on it was applied to several par- 
tial differential equations with variable coefficients containing some equations 
arising in the mathematical cosmology. More precisely, the new integral op- 
erator is suggested which transforms the family of the fundamental solutions 
of the Cauchy problem for the equation with the constant coefficients to the 
fundamental solutions for the operators with variable coefficients. The kernel 
of that transformation contains Gauss's hypergeometric function. 

This method was used in [51 to investigate in the unified way 

several equations such as the linear and semilinear Tricomi and Tricomi-type 
equations, Gellerstedt equation, the wave equation in Einstein-de Sitter space- 
time, the wave and the Klein-Gordon equations in the de Sitter and anti-de Sit- 
ter spacetimes. The listed equations play important role in the gas dynamics. 



elementary particle physics, quantura field theory in the curved spaces, and 
cosmology. For all above mentioned equations, we have obtained among other 
things, fundamental solutions, the representation formulas for the initial-value 
problem, the Lp — L^-estimates, local and global solutions for the semilinear 
equations, blow up phenomena, self-similar solutions and number of other re- 
sults. 

The starting point of our approach is the Duhamel's principle, which we 
revise in order to prepare the ground for generalizations. It is well-known that 
the solution of the Cauchy problem for the string equation with the source term 

/, /(a;,0eC-(M2), 

utt - ^ f{x,t) inR^, u(a;,io)=0, ut{x,to)=0 inK, (1.1) 
can be written as an integral 

/■* 

u{x,t) = / v{x,t;T)dT 

J to 

of the family of the solutions v{x, t; t) of the problem without the source term, 
but with the second initial datum 

Vtt-Vxx=0 inM^, v{x,t;t)^0, Vt{x,T;T) ^ f{x,T) in M . 

Our first observation is that we obtain the following representation of the 
solution of (|1.1[) 

u{x,t) = / dr w{x, z;t) dz , (1-2) 
Jto Jo 

if we denote ^ 

w{x- 1; t) - [f{x + i, t) + f{x - t, r)] , 

where the function w = w{x; t; r) is the solution of the problem 

Wtt-Wxx^O in M^, w{x,0;t) f{x,T), wt(x, 0;r)=0 in M . (1.3) 

This formula allows us to solve problem with the source term if we solve the 
problem for the same equation without source term but with the first initial 
datum. We claim that the formula (|1.2p can be used also for the wave equation 
with X G R", for all n e N. (See, e.g, [11].) More precisely, it holds also for the 
problem 

utt- f{x,t) inR"+\ u{x,to) = 0, ut{x,to) ^ in R" , 

with the function w ~ w{x; t] r) solving 

wtt - Aw = in R"+\ w{x,0;t) = f{x,T), wt{x,0;T) = in R". (1.4) 

Note that in the last problem the initial time i = is frozen, while in the 
Duhamel's principle it is varying with the parameter r. 
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The second observation is that in (|1.2p the upper hmit t — t of the inner 
integral is generated by the propagation phenomena with the speed which is 
equals to one. In fact, that is a distance function between the points at time t 
and r. 

Our third observation is that the solution operator G : f i — > u can be 
regarded as a composition of two operators. The first one 

WE : / I — >w 

is a Fourier Integral Operator (FIO), which is a solution operator of the Cauchy 
problem with the first initial datum for wave equation in the Minkowski space- 
time. The second operator 

K- : w I — > u 

is the integral operator given by (11.21) . We regard the variable z in (|1.2p as a 
"subsidiary time" . Thus, G = /C o WE and we arrive at the diagram: 



W£ 




« << 

u 

Our aim is that based on this diagram to generate a class of operators for 
which we will obtain explicit representation formulas for the solutions. That 
means also that we will have representations for the fundamental solutions of the 
partial differential operator. In fact, this diagram brings into a single hierarchy 
several different partial differential operators. Indeed, if we take into account the 
propagation cone by introducing the distance function <p(t), and if we provide 
the integral operator with the kernel K{t; r, b) as follows: 

K.[w\{x,t) ^2 j db j K{t;r,b)w{x,r;b)dr, x£W\ t>tQ, (1.5) 

Jto Jo 

then we actually can generate new representations for the solutions of the differ- 
ent well-known equations. Below we illustrate the suggested scheme by several 
examples. 

l" Klein-Gordon equation in the Minkowski spacetime. If we choose the kernel 
K{t; r, b) as 

Kit;r,b) = Jo{Vit-br-r') , (1-6) 

where Jo{z) is the Bessel function of the first kind, and if we choose the distance 
function as (j){t) = t, then we can prove (see Theorem 11.11 below) that the 
function 



u{x,t) = db J Jo (^V(t - 6)2 - r^^ w{x,r;b)dr, x e K, t > to, 
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solves the problem for the Klein-Gordon equation with a positive mass equals 
to 1 in the one-dimensional Minkowski spacetime, 

utt — + u ^ f{x,t) in u{x,to) — 0, ut{x,to)—0 inM, 

provided that w{x, r; b) is a corresponding solution of the problem for the 
wave equation in the Minkowski spacetime. We emphasis that the function 
w = 'w{x,t]b), with 6 regarded as a parameter, and the function u — u{x,i) 
solve different equations. This is a fundamental distinction from the Duhamel's 
principle. 

Now if we choose the kernel K{t; r, b) as 

K{t-r,b) = lo [^/{t - 6)2 - r^) , (1.7) 

where Io{z) is the modified Bessel function of the first kind, and the distance 
function as (/)(t) ~ t, then the function 

u{x,t) = db J /o w{x,r;b)dr, x G M", t > to, 

solves the problem for the Klein-Gordon equation with an imaginary mass in 
the one-dimensional Minkowski spacetime, 

utt - Uxx - u ^ f{x,t) inR2, u{x,to) = 0, ut{x,to)=0 inM, 

provided that w{x,r; b) is a corresponding solution of the problem ()1.3|) for the 
wave equation in the one-dimensional Minkowski spacetime. According to the 
next theorem the representation formulas valid also for the higher dimensional 
equations. The proof is by straightforward substitution. 

Theorem 1.1 The functions u — W7je(x,i) and uxrn{x,t) defined by 

une{x,t) = db Jo (rn^/{t - bf - r^^ w{x, r; b)dr, x G M", (1.8) 

ui,n{x,t) = db lo (771 7(1^^6)2^^) w{x,r;b)dr, x G K", (1.9) 

m = \M\, are solutions of the problems 

utt - Au + M'^u = f{x, t) in R"+\ u{x, to) = 0, ut(x, to) = in M" , 

with Af2 > and M'^ < 0, respectively. Here w{x,t;b) is a solution of 

Definition 1.2 The integral operator lll.5\) is said to be a generator of the so- 
lution operator for some equation if the operator G = K. o WE gives a solution 
operator for that equation. 
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2" Tricomi-type equations. The first example linkings to the operator with the 
variable coefficient is generated by the kernel K{t; r, 6) = E{0, t; r, b), where the 
function E{x, t; r, b) [TI] is defined by 



Eix,t;r,b) Cfc + - (x - r)^) ^ 

V'^'''m)+^{b)r-ix-r)^j ' ^'■'''> 

with 7 Cfe = (fc + i)-fc/(fc+i)2-i/('=+i), 2k = I e N, := [xp, and the 

distance function (p — (f>{t) is 

0W = fc^*'"-'. (1-11) 

while -F(a, &; c; C) is the Gauss's hypergeometric function. Here we assume that 
2fc € N but later on we consider the case of / G M. It is proved in p4 that for 
an integer non-negative I, for the smooth function / = f(x,t), the function 

u{x,t) = 2ci j dbj ((0(t)+0(&))2-r2) ^ 

solves the Tricomi-type equation 

utt-t^Au^ f{x,t) in Wl+^ :^ {{x,t)\x eK",t>0}, (1.12) 

and takes vanishing initial values 

u{x,0) = 0, ut{x,0)=0 inR". (1.13) 

3° The wave equation in the Robertson- Walker spacetimc: De Sitter spacetime. 
The next interesting example we obtain if we set K{t; r, b) — E{0, t; r, 6), where 
the function E{x, t; r, b) [13] is defined by 



E{x,t;r,b) := {{e-'' + e-'f - {x - rf) 

xFrii-l-(^l^^:!l!^(^) (114) 
^^2'2' '(e-* + e-^)2-(x-r)2 ' ' ^^''^^ 



and (/)(t) :— e *. It is proved in [14 that defined by the integral transform (II. 5p 
with the kernel (|1.14p the function 



u{x,t) = 2/ db I {{e-'' + e-*f - r^) ' 

Jq Jo 

fl 1 (e"*-e-^)2-r2\ , 
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solves the wave equation in the Robertson- Walker spaces arising in the de Sitter 
model of the universe (see, e.g. [8]), 

Utt-e-^'Au^ fix,t) mRl+\ 

and takes vanishing initial data (|1.13p . 

4P The wave equation in the Robertson- Walker spacetime: anti-de Sitter space- 
time. The third example we obtain if we set K{t; r, h) — E{0, t; r, b), where the 
function E{x,t;r,b) is defined by (see |15) ) 

E{x,t;r,b) := {{e'' + e'f - {x - rf)'^ 

while the distance function is (j){t) :— e* . In that case the function u = u{x,t) 
produced by the integral transform (jl.Sp with — Q and the kernel ()1.15p . solves 
the wave equation in the Robertson- Walker space arising in the anti-de Sitter 
model of the universe (see, e.g. ISj), 

utt - e^* Au = /(x, t) in M'|+\ 

Moreover, it takes vanishing initial values (I1.13p . 

The wave equation in the Einstein-de Sitter spacetime. If we allow negative 
; e M in (|LT0)) and, in that way, simplify the Gauss's hypergeometric function 
of the kernel of the integral transform, then we obtain another way to get 
new operators of the above described hierarchy. In fact, in the hierarchy of the 
hypergeometric functions F (a, &; c; C) there are functions which are polynomials. 
This is a case, in particular, of the parameter a — —m, where m € N. More 
precisely, if m G N, then I = ~4m/ {2m + 1) > —2 and F(— to, — m;l;C) = 
^(m-l)— (m+i-w) ^ ^„ that casc wc choose the distance function 
4>(t) — (2to -I- 1)^2;^ and the kernel K(t] r, b) as follows 

' m{m — 1) • • • (to + 1 — n) 



K{t; r, b) 



n=0 



(2m + l)2(t5T^ + 62^^)2 _ 

/ 11 \ ^ 

X ((2m + l)2(t27^ - 627^)2 - r^j . (1.16) 

Thus the integral transform JC allows us to write the representation for the 
solution of the equation 



utt-t 2™+iAu = / 



Moreover, in the hierarchy of the hypergeometric functions the simplest non- 
constant function is F(— 1,— 1;1;C) = 1 + C- The exponent / leading to the 
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function F (—1, —1; 1; () is exactly the exponent / = —4/3 of the wave equation 
(and of the metric tensor) in the Einstein & de Sitter spacetime. In that case 
of TO = 1 the kernel K(t; r, h) of (fTTB)) is K{t; r,b) = (9^2/3 + 952/3 _ j.2y 
Consequently, the function 

u{x,t)^ db —({3t^^^f + {3b^/^f-r^w{x,r;b)dr, (1.17) 

Jo Jo 18 V / 

a; G M", t > 0, solves (see [S]) the equation 

utt-t-^/^'Au^ f in Wl+\ (1.18) 

and takes vanishing initial data (jl.l3p provided that the function w is the image 
of /, that isw — yV£'(/). Because of the singularity in the coefficient of equation 
(|1.18l) . the Cauchy problem is not well-posed. In order to obtain a well-posed 
problem the initial conditions must be modified to the weighted initial value 
conditions. 

In fact, the operator of equation (jl.lSp coincides with the principal part 
of the wave equation in the Einstein & de Sitter spacetime. We remind that 
the Einstein & de Sitter model (EdeS model) of the universe was first proposed 
jointly by Einstein and de Sitter in 1932. It is the simplest non-empty expanding 
model with the line-element ds"^ = —dt'^ + a'^t^^^ [dx^ + dy^ + dz^). The covari- 
ant linear wave equation with the source term / written in these coordinates 
is 

The last equation belongs to the family of the non-Fuchsian partial differential 
equations. There is very advanced theory of such equations (see, e.g., [71 [TU]). 
but according to our knowledge the weighted initial value problem suggested in 
[5] (see (fn9l) below) is the original one. Assume that f{x,t) S C°°(M" x (0, cx))), 
and that with some £ > one has 

\d'^f{x, t)\ + \tdtd^f{x, t)\ < Cat'-^ for all x e M", and for small t > 0, 

and for every a, \a\ < [{n + l)/2]. It is proved in [5] that the function 

^P{x,t) = — db b(9t^/^ + 9b^^^ -r^)w{x,r;b)dr 

lot ./n .In 



solves the problem 



J tjjtt - f-4/3 A V + 21-^ A = f{x, t), i > 0, a; e R", 
] lim tip(x, t) = 0, lim {tipAx, t) -f 7/;(a;, t)) = 0, a; G E", 

provided that the function w is the image of /, that is w = WE{f ). 



(1.19) 



7 



2 The fundamental solutions of the operators 



In this section wc apply the method from the previous section to construct the 
fundamental solutions of the operators for the above listed equations. They are 
hyperbolic equations and therefore they have the fundamental solutions with 
the support in the forward or backward light cones. First we consider the string 
eqTiation. In order to write the fundamental solution with the support in the 
forward light cone we look for £ E such that 

£tt - £xx = 5{x - xo)S{t - to) in M^, suppf C {{x,t) \t>to, xeR} 

Then £{x, t; Xo,to) = £{x — xo,t — to; 0, 0) and if we denote by D{xo, to) the for- 
ward light cone D{xo, to) := {{x, t) S | |a; — xq] < {t — to)}- It is well known 
that £{x,t;xo,to) = if {x,t) € D{xo;to) and £{x,t;xo,to) = otherwise. We 
follows approach of the previous section and rewrite this fundamental solution 
in the following way 

£{x,t;xo,tQ) = H{t-to) ( ° £^*™B(x-xo,z)dz, 

Jo 

where H{t — to) is the Heaviside step function. The distribution 

f = i {5{x + t) + 5{x - t)} 

is the fundamental solution of the Cauchy problem for the string equation: 

^string _ ^string ^ (x, 0) = S{x), ""»(a;, 0) = . 

The string equation is partially, in direction of time, hypoelliptic, that implies 
^string g C'^ (R^-j)' {R^)). Hence, for every test function e C°°(M"), we have 

<£{x,t;-,to),ip{-)> = H{t-to) f '° <£''"'''ix-;z),^{-)> dz. 

Jo 

Thus, the forward fundamental solution of the operator is given by the integral 
transform of the fundamental solution of the Cauchy problem for the wave 
equation corresponding to the first datum. 

We can generate a class of equations which allows explicit representation 
formulas for the fundamental solutions. Indeed, if we provide the integral trans- 
form with a kernel as follows: 

.|0(t)-<^(to)| 

£{x, t; xo,to) = H{t - to) / K{t; r, 6)f "'"''^(x - xo,r)dr , 

Jo 

a; e M", t gR, then we get the representations for the fundamental solutions 
of the wide class of partial differential equations equations. In particular, if we 
plug in the integral transform the kernels used in the previous examples, then we 
obtain the corresponding fundamental solutions with the support in the forward 
light cone. 



8 



1° Klein-Gordon equation in the Minkowski spacetime. If we choose the kernel 
K{t; r, b) (|1.6I) and choose the distance function as (f>{t) — t, then it can be easily 
verified (see Theorem 12.11 below) that the distribution 

£{x, t- xo,to) = H{t - to) ^ Jo [^J{t-bY-r^) E^'^^^x - xo,r)dr , 

X e M", t E R., is the forward fundamental solutions for the Klein-Gordon 
operator with a positive mass equals to 1 in the Minkowski spacetime, 

{df - A + l) £{x,t;xo,to) ^ S{x - xo)S{t - to) in R"+\ suppf Q D{xo,to). 

provided that £^°-'"^{x,t) is the fundamental solution of the Cauchy problem 
corresponding to the first datum with the support at the origin, for the wave 
equation in the Minkowski spacetime. We emphasis that the distributions 
£{x,t]Xo,to) and £^'^^'^{x,t) solve different equations. 

If we now choose the kernel K{t; r, b) (|1.7p and the distance function as 
(j){t) — t, then the distribution 

£ix, t; xo, to) = H{t - to) ^ " lo [V{t^W^) f - xo, r)dr , 

X G K", t € R, is the forward fundamental solutions for the Klein-Gordon 
operator with an imaginary mass in the Minkowski spacetime, 

(a^ - A - l)£{x,t;xo,to) = 5{x - xo)S{t - to) inM"+\supp£ C D{xo,to). 

The following theorem can be easily proved by direct substitution. 

Theorem 2.1 The distributions £ne{x,t; xo, to) and £xm{x,t; xo,to) defined by 

£neix,t;xo,to) ^Hit-to) " Jo [m^{t - b)^ - r^') £'""'^(2: - xo,r)dr, 

£im{x,t;xo,to) ^ H{t~to) " lo [mVit - 5)2 - £^--^x - xo,r)dr, 

X G K", t G M, are forward fundamental solutions for the Klein- Gordon opera- 
tors with a real and an imaginary mass 

df~A + M^ inM"+\ 

with M2 > and M"^ < 0, respectively. Here m = \M\ > and £'^°-'"''{x,t) 
is the fundamental solution of the Cauchy problem corresponding to the first 
datum with the support at the origin, for the wave equation in the Minkowski 
spacetime. 

2° Tricomi-type equations. If we now choose the kernel K{t; r, b) (ll.lOp and the 
distance function as (jl.lip . then it is proved in jllj that £(x,t] xo,to) is the 
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forward fundamental solution for the Tricomi-type equation (I1.12p : 



£{x,t-XQM = 2ciH{t-to)l {{4>{t) + cf>{to)f - r^) 

a: e M", to > 0, with the support in the forward light cone 

D{xo,to) {{x, t) e M"+i \\x-xo\< (t){t) - <p{to)} • 

3° Klein-Gordon equations in the Robertson- Walker spacetime. The integral 
transform and, in particular, its kernel and the Gauss's hypergeometric function, 
open a way to establish a bridge between the wave equation (massless equation) 
and the Klein-Gordon equation (massive equation) in the curved spacetime. 
Indeed, if we allow the parameter 7 of the function F(7, 7; 1; z) to be a complex 
number, 7 £ C, then this continuation into the complex plane produces the 
fundamental solutions i5+(a:, xq, to) for the Klein-Gordon operator in the de 
Sitter spacetime as follows 

£+ix,t;xo,to) = 2H{t-to) r ' (ie-'^-y {{e-'° + e'')^ ~ tY''^'^' 

Jo 

where the distribution £^"-'"^{x,t) is the fundamental solution of the Cauchy 
problem for the wave equation, while the non-negative curved mass M > is 
defined as follows: ■= ^ — ^ 0- The parameter m is mass of particle. 
The fundamental solution f _ {x, t; xq, to) with the support in the backward light 
cone admits a similar representation. The fundamental solutions £+{x, t; XQ,to) 
and £-{x,t;xo,to) are constructed in [TJ for the case of the large masses m > 
n/2. The integral makes sense in the topology of the space of distributions. 
The fundamental solutions for the Klein-Gordon operator in the anti-de Sitter 
spacetime can be obtained by time inversion, t — ^ —t, from the fundamental 
solutions for the Klein-Gordon operator in the de Sitter spacetime. 

Moreover, the analytic continuation of this distribution in parameter M into 
C allows us to use it also in the case of small mass < m < n/2. Corresponding 
equation 

utt - e~2* A u - M^u = 0, 

can be regarded as Klein- Gordon equation with an imaginary mass. Equations 
with imaginary mass appear in several physical models such as c/)^ field model, 
tachion (super-light) fields, Landau-Ginzburg-Higgs equation and others. 

More precisely, for small mass < tti < n/2 we define the distribution 
£+ix,t;xo,tQ) by 

£+{x,t;xo,to) 
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= 2H{t-to)J (4e-*«-*)-*^((e-*« +e-*)2-r2) 

4" The above listed examples hint at some necessary condition on the pair 0, 
K{t; r, b) in order for that pair to produce a generator of the solution operator 
for some partial differential equation. 

Theorem 2.2 Assume that the integral transform m.5\) with the distance func- 
tion (j) and the kernel K{t; r, b) generate fundamental solution 

£{x, t; Xo,to) - H{t - to) / K{t; r, - Xo,r)dr , 



"'0 

X S M", i e M, of the partial differential equation 

n n 

Utt - ^ {x + ^(6,(a;, t)u)^^ + b{t)ut + c{t)u = / (2.1) 

ij — 1 2—1 

with the real- analytic coefficients. Denote Vi = Vi{t) and V2 = V2{t) two linearly 
independent solutions of the ordinary differential equation 

V" + b{t)V' + c{t)V = 0, Vi{Q) = 1 = F2'(0), F/(0) = = V2{0) . 

Then the function K{t; r, b) satisfies identity 

rm)-m\ vmv2{t)-Vi{t)V2{b) 

for all t > b> 0. 

Proof. For every function / G C°°(K x [0,oo)), which for any given instant 
i > has a compact support in x, the function 

nt M(t)-m\ 
v{x,t) — db K{t;r,b)w{x,r;b)dr 
Jo Jo 

where 'w{x,r;b) —< £^°-'"'^{x— ■ ,r), f{- ,b) >, solves the equation (|2.ip and takes 
vanishing initial data. It follows 

r ft Mit)-4'{b)\ / r \ 

V{t) := J v{x,t)dx^ J '^^ J K{t;r,b)ij w{x,r;b)dx j dr . 

On the other hand, 

w{x,r;b)dx ^ F{b), F{b) :^ [ f{x,b)dx. 
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Hence, 



rt ( .|0(t)-0(&)| \ 

V(t) = I F{h) / Kit- r, h) dr\dh. 



At the same time from the equation (|2.ip we obtain 

^V{t) + h{t)jV{t)+c{t)V{t)=F{t). 

Hence 

Thus, for the arbitrary function / G C°°(R x [0, oo)) for aU t one has 
\LJ^'' )\.L J^[t,r,0)ar vAb)V'(b)-V'(b)V,(b)) 



db^O. 



V,{b)Vi{b)-V({b)V2{b), 
The theorem is proven. □ 

Corollary 11) For the Tricomi-type equation and for the wave equations in 
the de Sitter and anti-de Sitter spacetime the following identities hold 

f — --k+r / / t^+i , 
t-b^2ck ^— + ^— - r2 



k 

















fc+i 


't'=+l I 






fc+i 











^ \ k+1 fc+1 ^ 



2 

r2 



xF -,- ^2 \dr, t>b>0, 



2' 2'^' (e* + e'')2-r2 
t > b, respectively. 

2) For the Klein-Gordon equation in the de Sitter spacetime we have b{t) = 0, 
c{t) = and the following identities hold. If the mass term M"^ is positive, 
then 

-| — e —k—i^j 

±smM{t-b) = 2j {Ae-'-y^' (^{e-" + e-'f ~ r^^ " 

/I 1 (e^'' — e^*)2 — r^x 
xF{- + iM, - + iM\l: ) — r 7775 Trldr. 

V2 '2 ' ' (e-^ + e-*)2-r2/ 
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If the mass term is negative, then 



/I 1 fe"'' - e~M^ - r^^ 

x^^(^HM|,2HM|;l;L___)__),.. 



These identities were used in |12) and |19| to prove a blow up phenomenon 
for the semilinear Tricomi-type equation and the Klein-Gordon equation in the 
de Sitter spacetime. 

3 Estimate of the tail inside of the hght cone. 
De Sitter spacetime 

In this section we consider the equation 

utt - e"^*Au = . (3.1) 

The forward and backward fundamental solutions for the operator of the last 
equation is constructed in [14) . By means of those fundamental solutions the 
fundamental solutions of the Cauchy problem are given as the Fourier integral 
operators in the domain of hyperbolicity, t > 0, with the data prescribed at 
t = 0, 

u{x,0)^Mx), ut(x,0) = (^i(x), .tgM". (3.2) 

The formula for the solution of this problem is given by Theorem 0.6 [14) with 
M = 0. More precisely, the solution of equation (|3.1I) with the initial data ipo, 
ifii e C^f (R"), prescribed at i = is: 

u{x,t) = e^v^g{x,(l}{t)) + 2 / v^g{x,(l}{t)s)Ko{(j){t)s,t)(j){t) ds 

Jo 

+ 2/ v^-^{x,(j){t)s)Ki{(j){t)s,t)(t){t)ds, a; € R", t > 0, (3.3) 

where (/)(t) := 1 — e^*. The function Vip{x,4>{t)s) coincides with the value 
v{x^ <j){t)s) of the solution v{x, t) of the Cauchy problem vu — Au = 0, v{x, 0) = 
^{x), Vt{x,Qi) = 0. The kernels KQ{z,t) and Ki{z,t) are defined by Ka{z,t) := 
— [^i?(z, i; 0, &)] and Ki{z, t) := E{z, t; 0, 0), respectively, where E{x, t; r, 6) 
is given by ([1.14p . Thus, 



Koiz,t) 

1 



[(l-e-*)2-z2]^(l + e-*)2 

, 9t 9. 1 / 1 1 (l-e-*)2-22 

+ l-e-^* + z^)-F ' 



t N /l 1 (l-e-*)2-z2 



2' 2' ' (1 



2 V 2' 2' ' (l + e-*)2 



,0<z< l-e-\ 
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and 

K,iz,t) = ((1 + - zr^F I; 1; , < . < 1 - e-\ 

It is important that the formula p. 31) regarded in the topology of the contin- 
uous functions of variable t with the values in the distributions space I?'(]R"), 
is applicable to the distributions (po,(pi G 2?'(R") as well. 

Recall that a wave equation is said to satisfy Huygens' principle if the solu- 
tion vanishes at all points which cannot be reached from the initial data by a 
null geodesic, that is there is no tail. An exemplar equation satisfying Huygens' 
principle is the wave equation in n-|- 1 dimensional Minkowski spacetime for odd 
n > 3. According to Hadamard's conjecture this is the only (modulo transfor- 
mations of coordinates and unknown function) huygensian linear second-order 
hyperbolic equation. Counterexamples to Hadamard's conjecture, which have 
been found, do not change the fact that Huygens' property is a very rare and 
unstable, with respect to the perturbations, phenomenon. It is natural, there- 
fore, to ask if there are other hyperbolic second-order equations which preserve 
Huygens' property approximately, in the sense that the tail which is left behind 
the wave front is comparatively small. For the equation in the de Sitter space- 
time Huygens' principle is not valid [14 . In this section we show the way how 
our approach can be applied to derive the pointwise estimates of the tail. 

We call the "tail" the part containing the integrals of p.3p . that is 

T{x,t) ■■= 2 f v^^{x,<j){t)s)Ko{<j){t)s,t)<l)[t)ds 
Jo 

+ 2 v^,{x,(j){t)s)Ki{(j){t)s,t)(j){t)ds, xeM", t>0. 
Jo 

Hence, T{x,t) = u{x,t) — eivtp„{x,<p{t)). The tail is of considerable interest in 
many aspects in the physics, and in particular, in the General Relativity [1], [5]. 

In this section we restrict ourselves to the case of one-dimensional x. For the 
one-dimensional wave equation in the Minkowski spacetime Huygens' principle 
is not valid, and, consequently, one can not anticipate it for one-dimensional 
equation in the de Sitter spacetime. But the last one reveals all difficulties and 
technical details allowing to overcome those difficulties in the case oi x € 
with n > 2. The results for the case of n > 2 will be published in a forthcoming 
paper. We start with simple example. 

Example. Let (Pq{x) — H{x), (pi{x) = 0, where H is the Heaviside step 

function. Thenu{x,t) = ^e^H{x+l-e-*)+H{x-l+e-^)]+J^~'' ' [H{x-z)+ 

H{x + z)] Kq{z, t) dz and, consequently, T{x, t) = H(x — z)Ko{z, t) dz + 

/g ^ H{x + z)Ko{z, t) dz. Consider the point xq such that < so < 1 — e~*. 
Then we have 

^(xo, t) = -62 -f / Ka(z,t)dz+ I Ko{z,t)dz^-+ / Ko(z,t)dz 
2 Jo Jo 2 Jq 
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while T{xo, t) = ^ — ^e^ + J^" Kq{z, t) dz = u{xo,t) — ^e^ . In particular, 



\T{xo,t)\ 



\u{xo,t) ~ T{xo,t)\ 



1 - e-*/2 - 2e-*/^ 



Ko{z,t) dz<2(l-e-*/2), 



lira 



lim 



\T{x,t)\ 



t^oo x^(l-e-t)- \u{x,t) - T{x,t)\ 



= 2. 



Thus, the tail dominates the huygensian part of the solution. If a;o > 1: then 
u{xo,t) = 1, while T{xo,t) = 1 - e*/^ 

Suppose now that initial data ipQ and ipi are the homogeneous functions. 



(3.4) 



<fo{x) = Co\x\-% ifiix) ^ Ci\x\-'' . 
The next theorem gives the pointwise estimate for the tail. 



Theorem 3.1 Consider the Cauchy problem for the equation LS.1\) . L3. ^|) with 
n — 1 and the data ^3.4^ , where a,b d (1/2,1). Then in the light cone emitted 
by the origin, that is on the set {{x,t) | < 1 — e~*, t > 0}, the solution, and, 
consequently, the tail, satisfy 



\T{x,t)\ = 



u{x,t) - icoe^ + 1 - e'X" + |x - 1 + e^^r" 



< |Co|C (1 + t) e^*e'^*(l + e*(l - \x\)Yl'^''' 
+ |Ci|C (1 + t) e-^*e^*(l + e*(l - \x\)fl^~'' . 

Proof. The representation of the solution of the Cauchy problem for the one- 
dimensional case (n = 1) of equation p.ip is given by Theorem 0.4[T3] with 
M = 0. More precisely, the solution u = u{x,t) of the Cauchy problem 

Utt - e~^*u^x ^ , u{x,0) ^ ipo{x) , ut{x,0) ^ ipi{x) , 

with ipo, (fi of p.4p . can be represented as follows 

1 



u{x, t) 



-62 



(^0(2: + 1-e ') + ipo{x-l + e 



[(po{x - z) + ipaix + z)\Ko{z,t) dz 



i-e-* 

ipi{x — z) + Lpi{x -\- z) Ki (z, t)dz . 
Then, it is evident that the solution of the Cauchy problem with Ci = is 

u{x,t) = Coe5[|2; + l-e-*r'^ + |x-l + e-V^ 



'\x - + \x + z\-'']Kn{z, t) dz . 
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Consider for x > the first term in the last integral It can be estimated as 
follows: 



\x-y\~''Ko{y,t)dy 



< 



\x-y\~''\Ko{y,t)\dy 



l-e-' 



\x-yr 



[(l-e-*)2-y2y(l + e-t)2_y2 



i)f{ 



1 1.^. (l-e-*)2-y^ 



2' 2' ' (l + e-*)2-j/2 
1 



1 1 (1 - - 



2' 2' ' (l + e-*)2-?;2 
If we denote z = e* and make a change y = e~*r in the last integral, then 

|x-y|-''|ifo(2/,i)|rfy 



|za; — r| " 



dy. 



[{z - 1)2 - r2]^(z + l)2-r2 



1 1 (z-l)2 



(z + l)2 -r2. 

1 1 (Z-1)2 -2, 



-1 + ^ j2^V"2'2'^'(z + l)2-r2j 



{z + 1)2 - r2 

Define two zones Zi(£,z) and Z2{e,z), as follows 

(z-l)2-r2 



dr . 



Zi{e,z) 
Z-,{s,z) 



(z,r) 
(z. r) 



„ < £, < r < z- 1 

(z + 1)2 - r2 - ' - - 

(z-l)2-r2 

£ < -, TT, ^, < r < z - 1 



(z + 1)2 - r2 

respectively. Then we split integral of (I3.5P into the sum /i + /2 where: 



(3.5) 

(3.6) 
(3.7) 



[(z - 1)2 - r2] ^(z + l)2-r2 



1 1 (z - 1)2 - 



2' 2' ' (z + l)2 



+ (z2-l + r2)iF(- 



1 1 (z-l)2-r^ 
2' 2'"^' (z + l)2 -r2 



dr. fc = 1,2. 



First, we restrict ourselves to the first zone Zi(£,z). We follow the arguments 
have been used in the proofs of Lemma 7.4 T4! and Proposition 10.2 [14 . In 
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the first zone we have 

/I 1 {z-lf-r^ 
V2'2' ' (z + l)2-r2 



1 1 [z-lf-r^ 

J- 1 



V 2' 2' 



(z + l)2-r2 



4 (z + 1)2 - r 

i-l (-^)^--% o 

4 (z + 1)2 - r2 



(Z-1)2 



(^ + 1) 



(z-l)2-r2 ^ 
+ 1)2 - r2 



We use the last formulas to estimate the term containing the hypergeometric 
functions: 



+ (.^-l + r^)-F(--,-;l;i-^^^) 



l(z-l)^ 



8(^+1)2 



-(r2 + 02_ i) + 2(-z2 + z)| 



(z- l)2-r^ 
,(z + l)2-r2 

In that zone the integral can be estimated by: 



(I- 



z-l 



h < C \zx-y\-'' 



1 



The last integral can be written as follows: 



1 + 



1 



(z + l)2-y2 



dy. 



\zx-y\-'' 



1 



^/z^/FTT^ 



1 + z- 



z + l-y 



dy 



= \z{l -x)- «r"|^(l + s)"'/' + + (3.8) 

Consider the first term of the last sum. Inside of the light cone |a;| < 1 — e~*, 
that is, \zx\ < z — 1, and one has 



1 

\/z Ji 



z(l-x)-sr"(l + s)-i/2 ds 

z(l-x) 



1 rz(i-x) 

(z(l-a;)-s)-«(l + s)-5rfs 

Vz Jl 

+ 4= r {s-z{l-x))-'^{l + srids. 
yz Jz(i-x) 



On the other hand 

r.z(l-rr) 



1 /-^(i-^) 1 

y=j^ (Z{l-X)-S)-%1+S)--^ 



^ ds= ^(l + z(l-a;))-" 
VZ 



17 



+ z(l - x)^ 



while 



^ / (s-z(l-a;))-"(l + s)-3ds 

Jz{l-x) 

-(l + z(l-a;))5-° 



2(1-2^) 

1 e''^^0Fr(l - a) 

7^ r(l-a) 



+ + "i^ -,a,-, 



1 + z 



2' l + 



Thus, 



^ iz(i~x)-sr« 



1 



: (is 



< 



z-5 (1 + z(i _ 2;))-'' -^^F (e*''^ - 1) (1 + z(l - x))^ 



-2V2F(-,a,-, ^ 

V2 2' l + z(l-a;) 



r(i-a) 
r(l-a) 



+2z-4(l + z)t(l + z(l-a;))-'^e""F I ^ a, ^ 



1 + z 



2' 1 + z(l - x) 



(3.9) 



Here the arguments of the hypergeometric functions satisfy the inequalities < 
i+z(i-x) ^ ^^'^ ^ ^ i+l{i-x) ■^(-'^ — a:) > 1, a: > 0, 2 > 1. Since 

- i - a) > we have 



-,a,-, 

'2' '2 l + 



< const for all z(l — a;) > 1 . 



Then, to estimate the last term of p.9p we use (4) of Sec.2.10 [2 : for all 

z{l - x) > 1, z > 1, X > 0, 



l + z 



'2' '2' l + z(l-a;) 



X 1 



l + z 



l + z 



1 + z(l - x) 

1-a 



l + z 



l + z{l -x) 



F l,-;2-a;l 



1 + z(l - x) 
1 + z(l - x) 



l + z 



where A, = ^2 = ^^^gy^^. Hence, 

^-|(l + ^)^(l + ^(l_^))-a^('l ;3 l + z 



2' ' 2' l + z(l-x) 
l + z 



1 + z(l - x) 
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+A2Z-^{1 + Z)5(l + z(l - x))-" 

1-a 



1 + Z 



1 + z(l - x) 



X 1- 



l + z 



1 + z{l- x) 



r 1,-2 — al 

' ' 2' ' 1 + z 



Since a > 1/2 one has 



1 l + z(l-x) 

l,-;2-a:l ^^-^ '- 

' ' 2' ' 1 + z 



< const 



1 + z 



1 + z(l - x) 



a-1/2 



for all z(l — x) > 1. It follows 



2' l + z(l- x) 



< C2-5(l + z(l_x))l/2-« 



Finally 



^ /"|z(l-x)-s|-'^^^ds < Cz-^(l + z(l-a;))i/2-°. 



From this inequality now we derive estimate for the second integral of 



|z(l ~x)~ L_ ds < Czi{l + z{l - , 



In the second zone we have 



(z-l)2-r2 
e < < 1 and 



< 



(z + 1)2 - r2 - (z - 1)2 - r-2 - e[{z + 1)2 - r^] 

and due to 15.3.10 of Ch.l5[I] we obtain 



■1 1 (z-l)2-r^ 



<C(l + lnz), foraU (z, r) e ^2(6, z). 



V2' 2' ' (z + l)2-r2 
This allows us to estimate the integral over the second zone: 

I2 < 



l(z,r)eZ2(e,z) 

< (1 + lnz) 



[{z - 1)2 - r2]^(z + l)2 -r^ 



(1 + In z) dr 



< (1 + lnz) 



''^ ((z + l)2-r2)3/2 

z-1 



dr 



z 



\zx — r\ 



dr 



((Z+1)2 -r2)3/2 

< (l + lnz)Vi/ \z{l-x)-s\-''{l + sy^/^ds 
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Then from the last inequahty and (|2.2I) we obtain, 

h < C(l + lnz)z^(l + z(l-a;))^/^~° 
Now consider the case of Cq = 0. Then 



T{x,t) = Ci / 
Jo 



l-e-* 



\x- z\~'' + \x + z\~'']Ki{z,t) dz, 



and 



\Tix,t)\ < Ci 



l-e-' 



[\x-z\-' + \x + z\-']{{l 
X F 



e-'?-z^)~^ 



1 1 (l-e-*)^-z^ 
,2' 2'"^' (l + e-*)2-z2 

Consider case of x > 0. Then make change of variable z = re^* 



dz . 



\Tix,t)\ < Cie'" / \e'x-r\-\ie'+l)^-r^) 



11 (e* - 1)^ - 
' 2'2' ' (e* + l)2-r2 



dr. 



Denote z = e*. Then 



\T{x,t)\ < Cz^ / \zx-r\-\{z + lf-r^) 



2' 2'^' (z + l)2-r2 



dr. 



Since 1/2 < 6 < 1, then for all z > 1 the following estimate 



z-l 



|za:-r|-''((z + l)^-r^r^i^ -,-;l; 



1 l.^_(z-l)2- 



2' 2' ' (z + l)2-r2 



dr 



< C(l + lnz)z"i/2(l + z(l-a; 



2-h 



is fulfilled. To prove the last estimate we rewrite the argument of the hyperge- 
ometric function as follows ^f^^'^]^yiZl-i = 1 - ^z+iyi^r^ . If r > ^(z + 1)2 - 8z, 
then i^^_^^yi_^.i > \ and < 1 - ^^^^^yi_^-i < \ for such r and z imply 



,11 4z 
F \ - -• 1- 1 

' 2'2' ' (z + l)2-r2 



< C. 



Hence we have 



< c 



\iz + lf-rT'^F[^y-;l;l- 
\zx-r\-\{z + lf -r^y^dr 



4z 



(z + l)2-r2 



dr 



< Cz"i/2 / |za; - rr''(z + 1 - r)-^dr 

Jo 

< cz-'/\i + z{i-x)y/^-'. 
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If r < ^{z + 1)2 - 8z and z > 6, then 8 < 8z < (z + 1)^ - < (z + 1)^, implies 

4z 



,11 4z 
F -,-;l;l 



2' 2' ' (z + l)2-r 



< C 



In 



(z + 1)2 - r2 



< C(l + lnz). 



Hence 



1 /I 1 4z 
|zx - r\-\[z + 1)2 - r2)-.F (^-, -; 1; 1 - ^ _ ) 



< I |za; -r|-''((z + 1)2 -r-')- 2 C(l + In z)dr 

rz-l 







< C(l + In z)z~^/2 / \zx-r\-''{z + l~ry^^dr 

Jo 

< C{1 + In z)z-i/2(i + _ a;))i/2-b _ 

Theorem is proven. □ 
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